Phonon-induced dephasing of chromium colour centres in diamond 
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We report on the coherence properties of single photons from chromium-based colour centres 
in diamond. We use field-correlation and spectral lineshape measurements to reveal the interplay 
between slow spectral wandering and fast dephasing mechanisms as a function of temperature. We 
show that the zero-phonon transition frequency and its linewidth follow a power-law dependence on 
temperature indicating that the dominant fast dephasing mechanisms for these centres are direct 
electron-phonon coupling and phonon-modulated Coulomb coupling to nearby impurities. Further, 
the observed reduction in the quantum yield for photon emission as a function of temperature is 
consistent with the opening of additional nonradiative channels through thermal activation to higher 
energy states predominantly and indicates a near-unity quantum efficiency at 4 K. 



PACS numbers: 61.72.J-, 81.05.ug 

Diamond plays a key role in a wide range of appli- 
cations in both electronics [Tj and photonics [2 today. 
Due to its wide bandgap it hosts a variety of optically 
active centres in the visible and the near infrared part 
of the spectrum [3J. The nitrogen- vacancy (NV) centre 
has attracted great attention due to its remarkable 
spin properties as a model system for quantum tech- 
nologies 0HS]. Other optically active centres such as 
silicon- vacancy [7\, nickel [8], xenon [9] and chromium 
[10) have also shown desirable photonic properties 
such as short lifetimes, predominant emission into the 
zero-phonon line (ZPL), and spectral tuning via DC 
Stark effect [TTl413| . However, the ZPL spectrum of 
diamond colour centres are typically broader than the 
radiative linewidth evidencing the influence of dephasing 
mechanisms on the optical transitions [T3]. These 
dephasing mechanisms can range from slow spectral 
wandering of the transition due to charge fluctuations 
in the environment to fast dephasing processes. The 
former can be remedied by direct feedback control on 
the transition, which has been demonstrated for NV 
centres [15) . The latter is irreversible and presents a 
fundamental limit to the single-photon coherence even 
with feedback. Therefore, it is essential to identify the 
source for these mechanisms, as well as the extent of 
their contributions to the spectral broadening of these 
centres. Here, we investigate the dephasing mechanisms 
influencing the chromium colour centres and report 
that the predominant sources for fast dephasing are 
direct electron-phonon coupling and phonon-modulated 
Coulomb coupling to nearby impurities, as manifested by 
the power-law temperature dependence of the transition 
frequency and linewidth. Further, the temperature 
dependence of the quantum yield reveals the existence of 
a thermal excitation mechanism to higher energy states 



similar to the neutral silicon-vacancy centres in diamond 

Chromium centres can be formed either by chemical 
vapour deposition (CVD) growth, where the etched sap- 
phire substrate acts as the source for chromium atoms 
|19) . or by implantation into bulk diamond [IT]. While 
the centres in bulk diamond exhibit narrower emission 
linewidth than those grown in nano- or microcrystals, 
they suffer from low photon collection efficiency due total 
internal reflection, as well as reduced photostability [T5] . 
We therefore concentrate on centres in microcrystals 
(size ^2 microns) for the investigation of dephasing 
mechanisms. The experimental setup used for the optical 
measurements is illustrated in Fig. 1(a). A temperature- 
controlled cryostat and a homemade confocal optical 
microscope are used to access the centres optically. The 
chromium centres are excited by a continuous wave laser 
at 700 nm and the emission (~ 780 nm) is collected into 
a single-mode fiber. A time-correlated photon counting 
setup is used for excited-state lifetime measurements, 
as well as intensity-correlation measurements. The 
excited-state lifetimes of chromium centres in bulk, 
micro- and nano-crystal diamonds can vary significantly 
[T71HE]- Figures 1(b) and 1(c) display exemplary results 
for these measurements performed at 4 K on a single 
chromium centre in a microcyrstal, and the value of 
1.7 ns extracted from the measurements is consistent 
with the range of values reported to-date. Based on this 
value, transform-limited photon emission is expected to 
result in a coherence time of 3.4 ns. A calibrated and 
scanablc Michclson interferometer and a spectrometer 
with 8-GHz resolution are used in parallel to quantify 
the degree of dephasing via the first-order coherence of 
the photon emission [see Fig. 1(a)]. 
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Figure 1: (Color online) (a) Schematic of the confocal micro- 
scope used for the experiments. BS is a beamsplitter and LP 
is a 740 nm longpass filter to remove residual laser reflection 
from luminescence by the chromium emitters, (b) Lifetime 
measured for the 773.3 nm transition line in figure la) inset. 
A value of 1.7 ns is extracted from the single exponential fit. 
(c) Autocorrelation measurement for the 773.3 nm line with 
g(2)(0) = 0.28. This value reflects the finite system response 
time and background emission from the crystal. The respec- 
tive fit again reveals an excited state lifetime of 1.7 ns. 



Figures 2(a) and 2(b) display the measured inter- 
ference Visibility V = [I max - I m in)/(I m ax + I m in) a s 

a function of relative time delay due to path-length 
difference for photons from two colour centres situated 
in two different microcrystals (labelled A and B, with 
ZPL at 756 nm and 770 nm, respectively). Filled black 
circles are the experimental data and the solid red 
curves are the theoretical fits. Common to both colour 
centres is the coherence timescale approximately 30 
times lower than the predicted transform limit. The 
coherence functions for these centres are, however, 
fundamentally different from each other. Centre A 
shown in Fig. 2(a) displays a predominantly Gaussian 
decay of coherence (the solid red curve is a Gaussian 
fit with T\i e — 62 ps), whereas for centre B shown in 
Fig. 2(b) the coherence decays exponentially (the solid 
red curve displays an exponential fit with T\/ e = 57 
ps). The exponential form of the coherence function 
for centre B indicates a dominant contribution from 
the irreversible dephasing mechanisms to the emission 
spectrum, while the Gaussian-like coherence function 
of centre A limits an exponential contribution to a 
timescale not shorter than 210 ps. This corresponds to 
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Figure 2: (Color online) First-order field autocorrelation mea- 
surement, g^ 1 ', performed at liquid helium temperature on 
centres A and B, as well as a chromium centre implanted in 
bulk diamond, (a) For centre A the visibility extracted from 
interference fringes at different time delays between the two 
arms of a Michelson interferometer (filled black circles) ex- 
hibits a Gaussian decay pattern with a T\i e = 62 ps (solid 
red curve). An exponential function is shown for compari- 
son (dotted black curve). The measured g^' for a chromium 
centre located in bulk diamond (filled blue circles) display a 
coherence time longer than 700 ps (solid red curve), (b) For 
centre B the extracted visibility (filled black circles) follows an 
exponential decay profile (solid red curve) with a coherence 
time of 57 ps. Again, a Gaussian profile (dotted black curve) 
is shown for comparison along with the g^> of the chromium 
centre in bulk diamond. 



750-MHz upper bound to the irreversible broadening 
for the linewidth for this centre. The relatively short 
coherence times exhibited by both emitters as well as 
the different profiles in spectra and coherence functions 
indicate that the immediate environment of the emitters 
has a varying impact on their photonic properties. In 
order to emphasize the effect of the environment on 
photonic coherence, similar measurements performed on 
a chromium centre implanted directly in bulk diamond 
with lower density of impurities is shown in both panels 
as filled blue circles. This centre shows a ZPL at 790 
nm. A coherence time longer than 700 ps is extracted 
for a Gaussian coherence function (solid red curve) - 
an order of magnitude longer than the coherence time 
measured for nanodiamonds. 
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The degree of optical dephasing can be revealed in 
first-order coherence as well as spectral measurements. 
If the line broadening is small compared to the natu- 
ral linewidth, the former is more accurate, whereas for 
strong dephasing, the latter provides a better means of 
determining emission linewidths. Figure 3(a) displays 
the characteristic broadening and red shift of the emis- 
sion spectrum of another chromium centre (labeled cen- 
tre C) as a function of lattice temperature between 9 K 
and 126 K. An analysis of the spectrum reveals that the 
Gaussian contribution remains constant for all tempera- 
tures, while the Lorentzian contribution is temperature 
dependent. This leads to a switch over from a predomi- 
nantly Gaussian-like lineshape at low temperatures to a 
Lorentzian lineshape at higher temperatures for the 10 
emitters investigated. This can be seen in the linear- 
log plots of the lineshape at two different temperatures 
displayed in the two insets of Fig. 3(a). The diamond 
bandgap itself will be modified as a function of temper- 
ature, following the Varshni law [20] [2T], and shallow 
centres will inherit this behaviour. Figure 3(b) displays 
the fractional change of the transition frequency (Aw/w) 
as a function of temperature, following a power law of 
the form Auj/uj oc T a with a=3.3. This frequency shift 
is distinctly stronger than the Varshni temperature de- 
pendence of bandgap energy for diamond, as is evident in 
Fig. [3jb) (black dotted curve). This reveals that all opti- 
cally relevant states are sufficiently far in energy from the 
band edges, as expected for deep defects in diamond. Fig. 
3(c) shows the temperature dependence of the linewidth 

7 of the Lorentzian contribution to the emission spectrum 
of centre C, following a power law of the form 7 oc 
with /3=3.7. These measurements were repeated for 10 
separate centres, and the average values of a and /3 are 
3.5 ± 0.3 and 3.4 ± 0.3, respectively. In the low temper- 
ature regime data suggests a temperature-independent 
dephasing mechanism is present for this centre limiting 
the linewidth to 4.1 GHz. This residual linewidth varies 
from centre to centre. We note however that a value of 
4.1 GHz is within the resolution of the spectrometer (~ 

8 GHz). 

A strong candidate for such temperature-dependent 
broadening mechanism is interaction with lattice 
phonons. The simplest model for electron-phonon cou- 
pling considers quasi-elastic scattering of phonons with 
a Debye density of states [22 from a non-degenerate 
transition of a centre. This model predicts a In- 
dependence for dephasing and a T 4 -dependence for the 
fractional change of the transition frequency. The strong 
deviation of the /3 values for chromium centres from 
the prediction by this model suggests the presence of 
an additional mechanism for temperatures below the 
Debye temperature Tr> An example of this is already 
seen in the NV linewidth, which exhibits a T 5 behaviour 
in ultrapure diamond due to Jahn- Teller effect in the 



excited states [23] • Further, experimentally determined 
values of a = 3.1 and /3 = 3 have been reported for 
the N3 colour centre related to nitrogen aggregates in 
diamond |24| . The particular temperature dependence 
was attributed to a pronounced deviation of the actual 
phonon density of states in diamond from the Debye 
model |25j . which was confirmed experimentally for 
the N3 colour centre using the vibrational structure in 
the N3 absorption spectrum. A numerical evaluation 
of linewidth and transition frequency shift using the 
modified phonon density of states in bulk diamond as a 
function of temperature according to Ref. [22] results in 
a = 3.1 and b = 3.7, which is in reasonable agreement 
with the values reported here. This model is plausible 
for chromium centres as well, although an experimental 
mapping of the phonon density of states is not feasible 
due to the small Huang-Rhys parameter [25] of these 
centres, which is less than 0.05 [18]. That said, the 
variation of a and (3 from centre to centre suggests that 
a mechanism governed more by local rather than global 
properties of the material may be playing a central role 
in the photon dephasing. 

An alternative model takes into account explicitly 
the impurity-rich nature of a material and incorporates 
the effect of phonons to the coupling of the centres 
(e.g. Coulombic and spin-dipolar) to nearby lattice 
impurities and other centres [27] , Consequently, a 
time-varying potential at the location of a centre leads 
to a T 3 dependence of the dephasing rate, valid for the 
temperature range T < T < Tjj, with Tjj the Debye 



temperature and To 



3/8 Tn. The coefficient c is the 



impurity concentration in the host crystal. Applying 
this general principle to the chromium emitters, an 
impurity concentration as high as 1 part-per-million 
gives To ~ 10 K and = 1850 K for diamond, so the 
temperature regime over which the chromium centres 
are investigated lies within this range. The presence of 
charges in the immediate environment of the emitter is 
already evident from the slow-wandering type decay of 
coherence exhibited by some of the centres (Gaussian 
coherence profile in Fig [2] (a)). This confirms the 
plausibility of such a dephasing mechanism. Further 
investigation of the temperature dependent dephasing 
as a function of concentration c of defects is necessary 
for an unambiguous identification of the dephasing type. 

Coupling of a colour centre to a vibronic mode also 
leads to reduced ZPL intensity with increasing tempera- 
ture [H[1I1[2H]: 



I (T) =/ o (0)exp 



-S'coth 



2k B T 
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2k B T / 

where Iq is the fraction of emission intensity decaying 
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Figure 3: (Color online) (a) Temperature dependence of an emission line associated with a chromium impurity for lattice 
temperatures between 8.6 K and 125.5 K, with temperature increased in steps of 15 K. The insets are spectra of the emitter 
plotted on a linear-log scale at 8.6 K and 65.2 K, respectively, and demonstrate the transition from a predominantly Gaussian 
to a predominantly Lorentzian lineshape. The red curves are Voigt fits and the black dotted curves are a Lorentzian curve 
at 8.6 K and a Gaussian curve at 65.2 K shown for comparison, (b) The shift of centre frequency follows a T 3 ' 3 power law 
(black solid curve). Varshni temperature dependence is shown for comparison (black dotted curve), (c) Transition linewidth 
as a function of temperature (blue circles). The linewidth follows a T 3 ' 7 power law (black solid fit), and a T 7 dependence is 
shown for comparison. 



into the ZPL, S is the Huang-Rhys factor, lo is the vi- 
bronic mode frequency, h is the reduced Planck constant, 
&b is the Boltzmann constant, and Jo is the zeroeth order 
Bessel function of the first kind. For chromium centres 
S < 0.05 [IB] and ?wj is on the order of about 10 meV, 
similar to other diamond centres [3J EH US]. Equation 
(1) applied to chromium centres yields 95% of the ZPL 
intensity to be sustained within the temperature range 
we study here. In strong contrast, Fig. 3(a) displays the 
significantly more rapid reduction of the ZPL intensity 
observed as a function of temperature. Such a strong 
reduction in intensity is uncharacteristic for S < 0.05 
and suggests an additional mechanism taking part. 

The quantum efficiency, 77, of a transition is determined 
by the presence of nonradiative channels connecting the 
excited and the ground states of a centre and is defined 
as 77 = P r /(P r + P nr ) where P r (P nr ) is the probability 
for a radiative (nonradiative) transition. Nonradiative 
channels contributing to P nr can become available as a 
function of temperature via thermal excitation from the 
excited state of the centre to higher energy nonradiative 
states. These states can either belong to the centre, or 



they could be formed by the impurity density in the vicin- 
ity of the emitters. Assuming that the selection rules do 
not affect the thcrmalization process, the quantum effi- 
ciency becomes temperature dependent: 

v(T) = - -^r- (2) 

1 + ge k B T 

The nonradiative decay process is assumed here to have 
a barrier with activation energy E t h, and the constant 
9 = Sth/ffex is the ratio of degeneracies between the ther- 
mally activated level and the excited state [l6j . The 
presence of a thermal excitation mechanism to additional 
states is further evidenced by the bunching behaviour 
in intensity autocorrelation measurements |10j , which 
appears predominantly for room-temperature measure- 
ments and is suppressed at low temperatures. Including 
this process, the overall intensity change of the ZPL as a 
function of temperature is given by [121 [Ml 130] : 
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Figure 4: (Color online) (a) and (b) ZPL intensity as a function of temperature for varying fitting parameters g and E t h, 
respectively, (c) and (d) Intensity of ZPL emission as a function of temperature for chromium emitters C and D, respectively. 
The side panels show the thermal coupling of the excited state to a reservoir of higher-lying states in the energy bandgap. We 
note that the centre in panel (c) is the same as that in Fig. 3. 



Fig. 4 (a) presents the calculated intensity of the ZPL 
according to equation 3 as a function of the number 
degeneracy of the state and the temperature (for a fixed 
value of the energy separation, E t h = 19.6). Similarly, 
Fig. 4 (b) presents the intensity as a function of E t h 
and temperature for a fixed value of the degeneracy, 
g = 19.0. Clearly, the intensity is sustained until the 
temperature is high enough to allow thermalization 
to be significant and the number of degeneracy deter- 
mines how rapidly the intensity is reduced beyond this 
temperature. Figures 4c and 4d present experimental 
results (blue circles) for two separate chromium centres 
superimposed on the Eq. 3 (solid black curves) using 
Eth and g as fitting parameters. The corresponding 
values for best fits are given in Table 1 along with those 
reported for the neutral silicon-vacancy centre |16| for 
comparison. The values for both g and E t h for centre C 
displayed in Fig. 4(c) are higher than those for centre 
D displayed in Fig. 4(d). The difference in energy 
separation and degeneracy of the thermally occupied 
states is not surprising given that the chromium emitters 
are known to be in an environment comprising a high 
density of impurity atoms such as additional nitrogen, 
oxygen and silicon [lOj [18]. While the model assumes 
thermal coupling to only one additional state, we 
interpret from the high values of g and E t h that the 
chromium emitters couple to an ensemble of states 
with an effective energy barrier E t h and an overall 
degeneracy ratio g. Intuitively, the number of states or 



the degeneracy in the ensemble is expected to increase 
with larger energy barriers -Eth, since the number of 
allowed energy states introduced by impurities around 
the chromium emitters should be higher closer to the 
valence band edge. This situation is depicted in panels 
(c) and (d) of Fig. 4. In principle, this model can 
be extended to include a quasi-continuum band with 
density of states g(E) formed by the impurities at a 
mean energy barrier E t h- Straightforwardly, P nr in 
Eq. (2) then becomes / g(E)exp(—^j i )dE. We note 
that using a Gaussian density of states of finite width 
around a mean energy barrier yields the same functional 
dependence, however independent measurements would 
be necessary to reveal the location of the energy barrier 
per emitter, supported by a density-functional theory 
calculations of the delocalized electronic states due to 
the impurities in diamond. The temperature-dependent 
model we use results in the values of 11% and 23% for 
the room temperature quantum efficiency for centre C 
and D, respectively. These efficiencies are consistent 
with previous reports on chromium centres at room tem- 
perature [17] and indicates that the strong suppression 
of thermal excitation allows quantum efficiencies at low 
temperatures, e.g. 4 K, of order unity. 

In summary, the temperature dependence of photon 
emission from single chromium centres in diamond 
revealed the interplay between slow and fast dephasing 
mechanisms and the nature of coupling to the impurities 
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g Eth S 

Centre C 19.0 ±8 19.6 ±3.5 meV 28.9 meV 0.05 
Centre D 4.3 ±2 4.4 ±1.7 meV 28.9 meV 0.05 
SiV° 1-3 5 meV 28 meV 1.5 



Table I: Fit parameters for the two chromium centres shown 
in figure and SiV centres for comparison. 

in the vicinity of the centres. Due to these coupling 
mechanisms, the photon coherence times remain about 
3.5 times in bulk diamond and 25 times in microcrys- 
tals below the theoretical value for transform-limit 
photon emission. Nonradiative losses due to thermal 
excitation to additional states are dominant at elevated 
temperatures, but the quantum efficiency of chromium 
centres approaches unity at low temperature operation. 
Deterministic implantation of these centres in ultrapure 
diamond, or their controlled incoporation in high purity 
CVD microdiamonds, is therefore needed to make these 
systems available as good-quality single photon sources. 
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